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Abstract. 

We present a numerical approach to investigate the relationship between magnetic fields and Faraday rotation 
effects in clusters of galaxies. We can infer the structure and strength of intra-cluster magnetic fields by comparing 
our simulations with the observed polarization properties of extended cluster radio sources such as radio galaxies 
and halos. We find the observations require a magnetic field which fiuctuates over a wide range of spatial scales 
(at least one order of magnitude) . If several polarized radio sources are located at different projected positions in 
a galaxy cluster, as is the case for A119, detailed Faraday rotation images allow us to constrain both the magnetic 
field strength and the slope of the power spectrum. Our results show that the standard analytic expressions 
applied in the literature overestimate the cluster magnetic field strengths by a factor of ~2. We investigate the 
possible effects of our models on beam depolarization of radio sources whose radiation traverses the magnetized 
intracluster medium. Finally, we point out that radio halos may provide important information about the spatial 
power spectrum of the magnetic field fluctuations on large scales. In particular, different values of the index of 
the power spectrum produce very different total intensity and polarization brightness distributions. 

Key words. Magnetic fields, Galaxies:clusters : general 



,!^ 1. Introduction 



The existence of cluster scale magnetic fields has been 
demonstrated in different ways. Some clusters of galaxies 
exhibit non-thermal synchrotron radio halos which extend 
up to Megaparsec scales. These diffuse radio sources con- 
stitute direct evidence for large scale magnetic fields as- 
sociated with the intracluster medium. Indirect evidence 
for cluster scale magnetic fields are provided by the obser- 
vation of Faraday rotation in radio sources located inside 
or behind galaxy clusters. For a review of these and other 
techniques see Carilli & Taylor (2002). 

A crucial issue concerns the effective strength and 
structure of cluster magnetic fields since different meth- 
ods of analysis give somewhat discrepant estimates for 
the field strengths. Due to their low surface brightness, 
radio halos have been studied so far with low spatial res- 
olution. This prevents a detailed observation of the small 
scale magnetic field geometry and intensity. However, us- 
ing minimum energy assumptions, it is possible to esti- 
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mate an equipartition magnetic field strength averaged 
over the entire halo volume. These estimates give equipar- 
tition magnetic field strengths of ~0.1 to 1 /iG (e.g. Feretti 
& Giovannini 1996, Govoni et al. 2001a, Bacchi et al. 
2003). In a few cases, clusters containing a radio halo 
show a hard X-ray excess emission. This emission could 
be interpreted in terms of Inverse Compton scattering of 
the cosmic microwave background photons with the rela- 
tivistic electrons responsible for the radio halo emission. 
In this case, the measures of the magnetic field strength 
(e.g. Fusco-Femiano et al. 1999, Rephaeli et al. 1999) in- 
ferred from the ratio of the radio to X-ray luminosities 
are consistent with the equipartition estimates. Because 
of the lack of spatial information of these X-ray observa- 
tions also in this case we obtain a field strength averaged 
over a large region. Indirect measurements of the magnetic 
field strength can also be determined in conjunction with 
X-ray observations of the hot gas, through the study of 
the Faraday Rotation Measure (RM) of radio sources lo- 
cated inside or behind clusters. By using simple analytical 
approaches, relatively strong magnetic fields, from ^5 /iG 
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up to the values of 30 //G have been found in coohng flow 
clusters (e.g., 3C 295, Perley & Taylor 1991 and Allen et 
al. 2001; Hydra A, Taylor & Perley 1993) where extremely 
high RMs have been revealed. On the other hand, signif- 
icant magnetic fields have also been detected in clusters 
without cooling flows: the RM measurements of polarized 
radio sources through the hot intracluster medium leads 
to a magnetic field of 2 — 8 /iG with patchy structures 
on small scales in the range 2 — 15 kpc (e.g. Feretti et 
al. 1995, Feretti et al. 1999a, Govoni et al. 2001b, Taylor 
et al. 2001, Eilek & Owen 2002). From a statistical study 
of several clusters, Clarke et al. (2001) obtained magnetic 
fields of 4 — 8 /xG, in agreement with the above findings. 

The magnetic field strength obtained by RM studies is 
therefore higher than the value derived from the radio halo 
data and from Inverse-Compton X-ray studies. However, 
as pointed out by Carilli & Taylor (2002 and references 
therein), all the aforementioned techniques are based on 
several assumptions. For example, the observed RMs have 
been interpreted until now in terms of simple analytical 
models which consider single-scale magnetic fields, while 
equipartition calculations in radio halos assume spatially 
uniform magnetic fields. 

Newman et al. (2002) demonstrated that the assump- 
tion of a single-scale magnetic field leads to an overestima- 
tion of the magnetic field strength calculated through RM 
studies. Moreover, detailed magneto-hydrodynamic cos- 
mological simulations (Dolag et al. 2002) suggested that 
the cluster magnetic fields may span a wide range of spa- 
tial scales with a strength that decreases with distance 
from the cluster center. 

Finally, It has been pointed out that in some cases 
a radio source could compress the gas and fields in the 
ICM to produce local RM enhancements, thus leading to 
overestimates of the derived ICM magnetic field strength 
(e.g. Feretti el al. 1995, Rudnick & Blundell 2003). No 
investigation in this respect has been attempted in this 
work. 

So far, very little attention has been given in the lit- 
erature to the determination of the power spectrum of 
the intra-cluster magnetic field fiuctuations. Only recently 
EnBlin & Vogt (2003) and Vogt & EnBlin (2003) by using a 
new semi-analytic technique showed that, for those cluster 
sources for which a very detailed RM image is available, 
the magnetic field power spectrum can be estimated. 

The aim of our work is to provide an alternative numer- 
ical approach for investigating the strength and structure 
of cluster magnetic fields through Monte Carlo simula- 
tions. The simulations results are compared to the data 
through an interactive technique. The algorithms are as- 
sembled into a software package named FARADAY. 

In particular, we will use FARADAY to study the 
power spectrum of the intra-cluster magnetic field from 
two different but complementary points of view. First, we 
show that both the mean and the dispersion of the RM 
can be used (and therefore should be used) to constrain 
not only the strength but also the power spectrum slope 
of the magnetic field fluctuations. In a consistent way, our 



simulations can be also employed to derive the expected 
de-polarization of cluster sources, such as radio galaxies 
and large scale halos, thus improving the estimates based 
on the RM analysis. 

Throughout this paper we will assume a simple power 
law form for the magnetic field power spectrum and com- 
pare the observed cluster radio source polarization and 
RM images with the simulated values over regions of com- 
parable size. The reader should be aware that many results 
presented in this work are infiuenced by the form assumed 
for the power spectrum and by the limitations of the nu- 
merical simulations. However, by considering the magnetic 
field strength as radial dependent and by modeling the 
magnetic field fluctuations within a multi-scale length ap- 
proach we made a new step towards more realistic cluster 
magnetic flelds models. As we will show, modeling of mag- 
netic fields as those we propose are necessary to interpret 
Faraday rotation effects in galaxy clusters and should be 
seen as a starting point for a new line of investigations. 

This paper is organized as follows: in Sect. 2 we give 
a brief description of the FARADAY tool. In Sect. 3 we 
discuss the three-dimensional multi-scale cluster magnetic 
fields models as implemented in the program and used 
in our simulations. In Sect. 4 we simulate the Rotation 
Measure. In Sect. 5 we simulate the beam depolarization 
of cluster radio galaxies and the morphology and depolar- 
ization of the radio halo emission. In Sect. 6 we present 
some comparisons of simulations with real data. Some dis- 
cussion and conclusions are presented in Sect. 7. 

In this work we adopt a Hubble constant 
Ho=50 km s^^Mpc^^ and a deceleration parameter 
qo = 0.5. 

2. Faraday rotation in clusters of galaxies and the 
FARADAY tool 

The polarized synchrotron radiation incoming from radio 
sources located inside or behind a galaxy cluster, experi- 
ences a rotation of the plane of polarization as it passes 
through the magnetized and ionized intra-cluster medium: 



*Ob.(!^) = *i„t + (c/iy)^ X RM 



(1) 



where ^obs(i^) is the observed polarization angle at a fre- 
quency v and '^int is the intrinsic polarization angle. 

The RM is related to the thermal electron density, Ue, 
and magnetic field along the linc-of-sight, i3[|, through the 
cluster by the equation: 



RM = 812 



/ 



UcB^^dl 



rad 



(2) 



where -B|| is measured in fiG, n^ in cm^'^ and L is the 
depth of the screen in kpc. 

The position angle of the polarization plane is an ob- 
servable quantity, therefore, RM images of radio sources 
can be constructed, by a linear fit to Eq. ^ The observed 
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degree of polarization intensity, Pobs(j^), can be signifi- 
cantly lower with respect to the intrinsic value, Pint, if 
differential Faraday rotation occurs within the observing 
beam. The observable data, RM and (de)polarization, to- 
gether with a model for the density distribution of the 
X-ray emitting hot gas can provide, in principle, impor- 
tant information on the cluster magnetic field strength and 
structure. However, there are a number of practical diffi- 
culties that complicate the applicability of this method. In 
particular, due to the random character of the cluster mag- 
netic fields, Eq.|21is not analytically solvable even for very 
simple distributions of the thermal electrons. Moreover, 
we will always observe a projection of the random field 
onto the two-dimensional sky, so we would need to de- 
project the power spectrum and the auto-correlation func- 
tion radio source RM images in order to reconstruct the 
magnetic field configuration. This may be a non-trivial 
process (Enfihn & Vogt 2003). 

For this purpose we developed the simulation tool 
FARADAY which allows an interactive, numerical ap- 
proach for investigating cluster magnetic fields. 

The thermal electron gas density of the screen is gen- 
erally not known, and we then assume a standard /3-model 
distribution: 



ne(r)=no(l-frVr,2)-W2 



(3) 



where r, no and Tc are the distance from the cluster center, 
the central electron density, and the cluster core radius, 
respectively. 

By following some simple steps it is possible to com- 
pute the effects of a cluster magnetic field model on the 
polarization properties of sources located inside or be- 
hind a cluster. Given a three-dimensional magnetic field 
model and the density distribution of the intra-cluster gas, 
FARADAY calculates the RM image by integrating Eq. |21 
numerically. 

The beam depolarization is expected to be particularly 
strong if the radio source is observed with a beam larger 
than the minimum magnetic field scale length. This effect 
can be studied with FARADAY as well. Starting from the 
intrinsic polarization images of Pj^t and ^int, FARADAY 
can simulate the expected beam depolarization as a func- 
tion of wavelength. One of the key feature of FARADAY is 
that it allows us to study the effect of magnetic field mod- 
els that cannot be investigated analytically. In this work 
we use FARADAY to investigate the effects of a turbulent 
and isotropic magnetic fields which fluctuate over a range 
of spatial scales. 

3. A multi-scale magnetic field model 

Detailed images of the polarized emission are now avail- 
able for a significant number of radio sources seen in 
galaxy clusters and indicate that the RM fluctuates down 
to linear scales as low as 10 kpc or less (e.g. Carilli & 
Taylor 2002 and references therein). If interpreted in terms 
of an external Faraday screen, the dispersion, ctrm, and 



mean, (RM) , of the RM fluctuations can be used to con- 
strain the cluster magnetic field strength. 

The effect of Faraday rotation from a tangled magnetic 
field has been analyzed by several authors (e.g. Lawler & 
Dennison 1982, Tribble 1991, Feretti et al. 1995 and Felten 
1996). In the simplest ideal case, the screen is made of 
cells of uniform size, electron density and magnetic field 
strength, but with a field orientation at random angles in 
each cell. The observed RM along any given line of sight 
is then generated by a random walk process involving a 
large number of cells of size Ac. The distribution of the 
RM is Gaussian with zero mean, and a variance given by: 



^2 



(RM^) = 812^Ac / {iieBnydl 



/' 



(4) 



As we will show in Sect l4.4l since the density profile 
of the ionized gas can be obtained by X-ray observations, 
the intergalactic magnetic field strength can be estimated 
by measuring ctrm from spatially resolved RM images of 
radio sources in the cluster if Ac is known. 

In many cases the observed RM distributions are 
nearly Gaussians, suggesting an isotropic distribution of 
the field component along the line-of-sight. However, many 
RM distributions show clear evidence for a non-zero 
(RM) if averaged over areas comparable with the radio 
source size which is typically 50 x 50 kpc^. This same 
property is observed for clusters located at high galactic 
latitudes and in general cannot be entirely attributed to 
the contribution from our own Milky Way Galaxy. These 
(RM) offsets are more likely due to fluctuations of the clus- 
ter magnetic flelds on scales greater than 50—100 kpc, i.e. 
considerably larger of those responsible for the RM disper- 
sion. The random magnetic field must therefore be both 
tangled on sufficiently small scales in order to produce the 
smallest structures observed in the RM images and also 
fluctuate on scales one, or even two, orders of magnitude 
larger. For this reason, it is necessary to consider cluster 
magnetic field models where both small and large scales 
structures coexist. To accomplish this we simulated in a 
cubical box a three-dimensional multi-scale magnetic field. 

Following the approach proposed by Tribble (1991), 
the field is constructed by selecting a power spectrum for 
the vector potential A and choosing the Fourier compo- 
nents A(k) accordingly. 

For all grid points in the wave-numbers space the am- 
plitude and phase, respectively A and 0, of each compo- 
nent of A(k) are randomly drawn from the distribution: 



P{A, (j))dAd(j) 



A 



27r MJ2 



exp 



A' 



2\Ak\y 



] dAd(j) (5) 



i.e. A is extracted from a Rayleigh distribution and cf) is 
uniformly distributed between and 27r. 
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We adopted a power-law power spectrum^, \Ak\ 



the form: 



lAkl^ock-'^. 



The magnetic field: 

B(k) =ikx A(k) 



of 



(6) 



(7) 



is then transformed back into the real space using a three- 
dimensional Fast Fourier Transform (Press et al. 1986). 
This automatically ensures a tangled, divergence free, 
multi-scale magnetic field model whose components have 
a power spectrum which follows: 



\B,f 



G„fc 



(8) 



where n = ^ — 2 while C^ is the power spectrum normal- 
ization. 

The power spectrum © represents the magnetic field 
energy (erg cm^) associated with each wave-number. 

The assumed form for the vector potential results in a 
magnetic field whose components, Bi, have a Gaussian dis- 
tribution with {Bi) = and dg, — {Bf). The distribution 

of B itself is therefore a Maxwellian with (B) = 2\ - ctr 



andfiB^^^ CTs.. 

From the Parseval's Theorem, j \B,\'^dV = j \Bk\^(fik, 
it follows that the magnetic field energy density averaged 
over the volume, V , is proportional to 



(B)2 ^ ^ . ln(fc^ax/fcmin) 

for n = 3, and 



ri2 h^-n _ jL3-n 
/D\2 „, ^n '^max '*'min 



(9) 



(10) 



for n 7^ 3. Where the wave- numbers k^^ 
correspond, respectively, to the minimum and maximum 
spatial scale of the magnetic field. 

3.1. A magnetic field strength with a radial decrease 

Important clues about the radial gradient of the magnetic 
field strength in clusters of galaxies are provided by obser- 
vations of clusters hosting a radio halo and by magneto- 
hydrodynamic cosmological simulations. The spatial cor- 
relation found in some clusters between the X-ray and the 
radio halo brightness (Govoni et al. 2001c) suggests that 
thermal and non-thermal (relativistic particles and mag- 
netic fields) energy densities could have the same radial 
scaling. In the radio halo of the Coma cluster (Giovannini 
et al. 1993) a steepening of the synchrotron spectrum is 

^ Note that throughout this paper the power spectra are 
expressed as vectorial forms in fc-space. The one-dimensional 
forms can be obtained by multiplying by 4Tvk^ and 2Tvk respec- 
tively the three and two-dimensional power spectra. 



observed from the cluster center outward. Such a spec- 
tral behavior is expected in some halo formation mod- 
els which consider a radial decrease of the cluster mag- 
netic field strength (e.g. Brunetti et al. 2002). Moreover, 
detailed magneto-hydrodynamic cosmological simulations 
(Dolag et al. 2002) lead to magnetic fields whose strength 
decreases with the distance from the cluster center. 

Therefore in this work we considered a magnetic field 
whose strength decreases from the cluster center according 
to: 



(B)(r) = (B)o-(l + rVre^)- 



(11) 



where (B)o is the mean magnetic field at the cluster cen- 
ter. Given the electron density profile in Eq. 13 we have 
(B) oc Tic . Therefore the case ji = (3/2 corresponds to 
a magnetic field whose energy density decreases from the 
cluster center as the gas energy density while ^ = 2fi/'i 
corresponds to a magnetic field frozen in the matter. 

By combining RM and X-ray data it is possible to 
obtain an estimate of the index fj, (Dolag et al. 2001). 
The X-ray surface brightness, S^, in a case of the /3-model 
(Eq. I^l, is proportional to 



^.oc(l + ri/r,2)-3^+4 



(12) 



while, by substituting Eq. ^^in Eq. ^ the RM is propor- 
tional to 



(RM2)i/2o,(i + r^/r2)-f(/3+M)+4 



(13) 



where r±_ is the projected distance from the cluster center. 
By comparing the two line-of-sight integrals, one finds that 
the index /i is related to the slope, a, of the correlation 
(RM^)^/^ oc S" and to the density index, /3, through 



^l= (2a - 1) • (/3 - 1/6) 



(14) 



Thus, for a constant magnetic field {^ — 0) the slope of 
the RM— S'x correlation should be a = 0.5 while a steeper 
slope would imply /i > 0. 

Formally, a decreasing magnetic field strength, as in 
Eq. ^2 should be achieved by convolving the vector po- 
tential components with the shaping profile in the Fourier 
space. A more practical, although less rigorous, way to 
obtain this, consists in shaping the magnetic field to the 
desired profile directly in real space. We checked that in 
our cases the two approaches give minimal differences and 
for simplicity we adopted the latter. It should be noted 
that the magnetic field power spectrum in both cases is 
given by the convolution of the power-law in Eq. [S] with 
the Fourier transform of the shaping function. However, 
since the core radius of the shaping function is a consider- 
able fraction of the simulated volume the convolution does 
not significantly alter the form of the power-law in Eq. |S1 

4. Simulated Rotation Measures 

Using the cluster magnetic field model described in Sect. 3, 
we studied the behavior of the RM for different values of 
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Table 1. Parameters adopted in the simulations. 



Grid size 
Cellsize 



512^ pixels 
1 pixel=3 kpc 



Core radius 
Central density 
Beta 



r(.=400 kpc 

no=10' 

/3=0.6 



no=10~^ cm"^ 



Central Mean magnetic field 
Radial profile slope 
Magnetic field minimum scale 
Magnetic field maximum scale 
Power spectrum spectral index 



{B)o = l ^^G 
^ = 0,0.3 

Amin=6 kpc 

Amax=768 kpc 
n=2, 3, 4 



the spectral index of the power spectrum in the range 
from n = 2 to n = A. The magnetic energy density scales 
as duB/dk oc fc-^jSfcp, thus n ~ 3 implies that the field 
is scale-invariant, i.e. the energy density per logarithmic 
wavenumber interval is constant. It follows that for n < 3 
the magnetic energy density is larger on the smaller scales 
while for n > 3 most of the magnetic energy density is on 
the larger scales. We now explore all these possibilities. 

4.1. Simulated RM images 

We simulated the expected RM images corresponding to 
the three different values of the power spectrum spectral 
index mentioned above in the case of a cluster with typical 
size and gas density parameters. The computational grid, 
gas density and magnetic field model parameters used for 
these specific simulations are listed in Tab. ^ while the 
RM images we obtained are shown in Fig. 1. Each RM 
image shows a field of view of about 1.5 x 1.5 Mpc while 
the cluster core radius (indicated by the circles in Fig. 1) 
is taken to be 400 kpc. We performed the integration of 
Eq. 13 from the cluster center up to 1.5 Mpc (i.e. 3.8 Tc) 
along the line-of-sight. To do this we located the cluster 
center in the farthest side of the computational grid. The 
three magnetic field power spectra are normalized to gen- 
erate the same average magnetic field energy density. The 
average field at the cluster center is (B}o ~ 1 ^G and 
its energy density decreases from the cluster center as the 
gas energy density profile, i.e. fi = 0.3 (see Eg. ITT)) . This 
corresponds to a mean magnetic field of (B)~ 0.46 fiG 
within a spherical volume with a radius of r = 3rc. 

Since the size of the computational grid is finite only 
a limited range of spatial scales can be investigated. The 
grid is large enough to contain a cluster of typical size 
while keeping the necessary resolution to sample reason- 
ably well the magnetic field structure on a spatial scale^ 
range from Amjn = 6 up to A^ax — 770 kpc. Here we 
adopted the maximum range of spatial scales permitted 
by the computational grid. Observations reveal RM fluc- 



^ Here we refer to the length A as to the magnetic field re- 
versal scale. In this way, A corresponds to a half-wavelength, 
i.e. A = 0.5-(27r/fc)- 



tuations on scales as low as 10 kpc supporting this choice 
for Aniin- The exact value of Amax is more uncertain, how- 
ever the presence in some cluster of galaxies of wide radio 
halos may indicate that at least in these cases the cluster 
magnetic field is spread over Mpc distances. Simulations 
with a different choice of of Amax are presented in Sect 14. 31 
The resulting two-dimensional power spectra of the 
RM images are shown in the top- left panel of Fig. 1. They 
have the same slope as their parent's three-dimensional 
magnetic field power spectra and they cover an equivalent 
range of spatial scales. However, it is evident from Fig. 1 
that the same cluster magnetic field energy density will 
generate different magnetic field configurations with cor- 
respondingly different RM structures, for different values 
of the power spectrum spectral index. 

4.2. RM radial profiles 

RM images of a whole cluster, such as those shown in 
Fig. 1, cannot be observed in real cases. However, it is 
relatively easy to measure the RM dispersion and mean 
(crm and (RM)) in limited regions by observing radio 
sources located at different projected distances from the 
center of a cluster. The ctrm and (RM) values found in 
radio sources, are usually calculated on their project sur- 
face which is typically about 50 x 50 kpc^ in size. Since we 
need to compare the simulation results with these observ- 
ables, it is necessary to calculate the simulated values of 
crm and (RM) in regions with the same projected area. 

Fig. 2 shows the profiles of ctrm (left column panels), 
|(RM)| (central column panels) and |(RM)|/crRM (right 
column panels) as a function of the projected distance 
from the cluster center, calculated in boxes of 50 x 50 
kpc^ in size. Due to the random nature of the simulated 
magnetic field, the values of ctrm and |(RM)| at a given 
radial distance vary from box to box. We calculated the 
mean and dispersion of ctrm and | (RM) | values from all 
the boxes located at the same projected distance from the 
cluster center. The lines and the shaded regions plotted 
in Fig. 2 represent the mean and dispersion of ctrm and 
(RM)|, respectively. 

The size of the sampling region used to trace the 
RM profiles, falls inside the power spectrum range scales. 
Therefore, ctrm is explained by the fact that the magnetic 
field fluctuates down to scales smaller than the size of the 
sampling regions. The most striking feature appearing in 
Fig. 2 is the significant amount of (RM) predicted by the 
simulations for any of the three spectral indexes adopted. 
This is due to fluctuations of the cluster magnetic field on 
scales comparable or greater than the size of the sampling 
regions. This is not an unexpected result, but the numer- 
ical simulations provide a quantitative evaluation of this 
effect. 

As expected, both ctrm and |(RM)| decrease from 
the cluster center outward, and the profiles correspond- 
ing to a constant magnetic field strength are the flattest. 
Moreover, the ratio of | (RM) | to ctrm as a function of dis- 
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Fig. 1. Simulated RM images for different values of the magnetic field power spectrum spectral index n. The three 
power spectra are normalized to have the same total magnetic field energy which is distributed over the range of 
spatial scales from 6 kpc up to 770 kpc. The average field at the cluster center is (B)o = 1 /iG and its energy density 
decreases from the cluster center according to B^ oc nc{r), where nc{r) is the gas density profile. Each RM image 
shows a field of view of about 1.5 x 1.5 Mpc while the cluster core radius (indicated by the circle) is 400 kpc. We 
performed the RM integration from the cluster center up to 3.8 core radii along the line-of-sight. The two-dimensional 
power spectra of the simulated RM images are shown in the top left panel. They have the same slope as their parent 
magnetic field power spectra and they span an equivalent range of spatial scales - see text for details. 
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Fig. 2. Radial profiles obtained from the RM simulations described in Sect. 01 Left, central and right panels show 
the radial trends of ctrm, KRM)! and |(RM)|/crR,M, respectively, for three values of the index of the power spectrum 
n. Solid lines correspond to an average field at the cluster center of (B)o = 1 /iG whose energy density decreases 
from the cluster center according to B^ ex nc(r), where nj(r) is the gas density profile. The shaded regions indicate 
the l-rms statistical fluctuations due to the random nature of the magnetic field. Dotted lines represent the profiles 
corresponding to a constant magnetic field strength of (B)o = 1 ^G all over the cluster. Values have been sampled in 
50 X 50 kpc^ regions. 



tance from the cluster center is constant, i.e. both quan- 
tities have the same average radial gradient. 

While both ctrm and |(RM)| increase linearly with the 
average cluster magnetic field strength, (B), and the cen- 
tral gas density (see Eq. |21), the ratio |(RM)|/(Trm de- 
pends only on the magnetic field power spectrum slope, 
n, for a given range of fluctuation scales. This means that 
the comparison of RM data of radio galaxies embedded in 
a cluster of galaxies with simulated profiles, can be used 
to infer the strength and the power spectrum slope of the 
cluster magnetic field. 

4.3. The dependence of RM on the magnetic field 
largest scale structure 

To trace RM profiles shown in Fig|21we adopted the largest 
magnetic field fluctuation scale permitted by the computa- 
tional grid size. In this section we consider different values 
for Amax and we calculate the corresponding variations of 
the simulated RM. 



We performed an additional series of simulations where 
we decrease the value of Amax from 768 kpc down to 6 kpc. 
We keep the minimum magnetic field scale Amin fixed at 6 
kpc and allowed A^ax to vary, since, given the limited size 
of the radio source RM images, the former can be better 
estimated from the data with respect to the latter. 

The variation of ctrm , | (RM) | and their ratio as a func- 
tion of Ainax is shown in left, middle and right panel of 
Fig. 3, respectively. 

Note that ctrm and | (RM) | in middle and left panels of 
Fig. 3 have been normalized to their values at Aniax=768 
kpc. Doing this we obtained a scaling factor for the profiles 
shown in Fig. 2 which depends only on the power spec- 
trum slope n and on Amax for a fixed Amin . We found that 
for n=3, and especially n=4, ctrm increases until Amax is 
larger than the size of the sampling region in which the 
statistic is calculated (i.e. 50 x 50 kpc^) while it stays 
almost constant for n=2. When Amax becomes smaller 
than the sampling region, ctrm starts to decrease. Instead, 
|(RM)| monotonically decrease by decreasing Amax and 
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has been fixed at Amin = 6 kpc while ctrm and | (RM) | values have been sampled in 50 x 50 kpc-^ regions. 



it goes to zero as A^ax -^ Amin for all the three spectral 
indexes. 

The right most panel of Fig. 3 shows the behavior of 
|(RM)|/(TptM as a function of A,„ax for the three different 
slopes of the magnetic field power spectrum. In general we 
find the ratio |(RM)|/ctrm decreases with A^ax- 

If Amax is significantly larger than about 50-100 kpc, 
a given value of |(RM)|/crRM can be explained by differ- 
ent combinations of parameters: the higher n, the lower 
Amax- This could be considered as a serious degeneracy at 
a first glance. However, it is important to bear in mind 
that the resulting RM images, albeit characterized by the 
same value of |(RM)|/(Trm, still have very different spa- 
tial power spectra and thus different correlation lengths. 
Therefore, if good RM images are available for the cluster 
radio sources, it is possible to break the degeneracy be- 
tween n and Amax, as will be shown for the case of A119 
in Sect.Hm 

As Amax becomes lower than about 50-100 kpc, the 
ratio |(RM)|/(7rm falls below 0.2-0.3 and its depen- 
dence on the power spectrum slope vanishes. In this case 
all the three spectral indexes predict almost the same 
|(RM)|/(Trm which can be then used to constrain Amax- 

In any case, observed values of |(RM)|/(7rm greater 
than unity can only be explained by steep magnetic field 
power spectra. For example, values of |(RM)|/crRM> 1 rule 
out n = 2 for any value of Amax- 

4A. Comparison with single-scale analytical magnetic 
field model 

The bulk of cluster magnetic field strength estimates ap- 
pearing in the literature have been calculated so far by 
using an analytical formulation based on the approxima- 
tion that the magnetic field is tangled on a single scale Ac. 
In this formulation, by considering a density distribution 
which follows a /3-profile, the following relation (hereafter 
analytical formula) for the RM dispersion with the pro- 



jected distance from the cluster center, r±, is obtained by 
integrating Eq. 4: 



ctrm(?'j_) = 



KBnorc' A 



1/2.1/2 



/r(3/3-0.5) 



(l+r2^/r2)(6/3-i)/4W r(3/3) 



(15) 



where F is the Gamma function. The constant K de- 
pends on the integration path over the gas density distri- 
bution: K = 624, if the source lies completely beyond the 
cluster, and K = 441 if the source is halfway the cluster. If 
the magnetic field scales with the distance from the cluster 
center as in Eq. ^2 P and B in Eq. El must be replaced 
by /3 -I- /i and Bq, respectively (see Dolag et al. 2001). 

It is of interest to compare the results of our simu- 
lations with the predictions of the analytical approach. 
However, since in this work we consider a magnetic field 
which fiuctuates over a range of scales extending at least 
one order of magnitude, it is necessary to discuss which is 
the value for Ac we should adopt in EalT^ 

Here we consider three possibilities: i) Ac is equal to 
the minimum field scale-length Amin, ii) Ac is equal to the 
field correlation length'^ scale Ab^ and iii) Ac is equal to 
the RM correlation length scale Arm. 

In Fig. 21 we compare the ctrm trend expected from the 
FARADAY simulations (solid line) with the predictions 
of the analytical formulation for the three scale lengths 



The magnetic field correlation length is defined as 



Ab, = 



Wz(0) 



(16) 



where ■u;z(r) is the spherically averaged autocorrelation func- 
tion of the magnetic field component along the line-of-sight. 
The RM correlation length is defined as 



Arm 



/o°°CRM(r)dr 
Crm(O) 



(17) 



where Crm('") is the radially averaged RM autocorrelation 
function. 
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(dashed and dotted lines). The simulated trend corre- 
sponds to a magnetic field power spectrum with n=2, 
Amin = 6 kpc and A„iax = 768 kpc. This model has a 
magnetic field correlation length of Ab^ — 16 kpc and a 
RM correlation length of Arm — 69 kpc. 

High-resolution RM maps of extended radio galaxies in 
clusters permit to detect patchy structures on small scales 
with typical size in the range 2-15 kpc (e.g. Feretti et al. 
1995, Feretti et al. 1999a, Govoni et al. 2001b, Taylor et 
al. 2001, Eilek & Owen 2002). The cluster magnetic fields 
strength reported in these works has been usually calcu- 
lated by considering Ac in Eq.^jto be equal to the small- 
est patchy structures detectable in the RM maps. This is 
the motivation for considering A^. = A,nin. In this case, by 
assuming the same magnetic field strength, the single-scale 
magnetic field approximation leads to a lower ctrm with 
respect to the prediction of our simulations. As a result 
the magnetic field strength inferred by the analytical for- 
mula is systematically overestimated by a factor of about 
2 with respect the magnetic field obtained by assuming a 
multi-scale magnetic field with a power spectral index of 
n = 2. 

On the other hand, since for a broad power spectrum 
Arm > Ab^ (see Vogt & Enfilin 2003), the assumption of 
the RM correlation length for Ac in Eq. El leads to under- 
estimate the magnetic field strength by the same factor. 

The best agreement between the simulated and ana- 
lytical ctrm trends is found for Ac — Ab^. In fact, this 
should be the proper scale length to use in Eq. E| 

To summarize, the analytical formulation leads to a re- 
liable estimate for the average magnetic field strength pro- 
vided that the magnetic field correlation length is known. 
Hoverer, since Ab^ is n dependent, we need to estimate 
the power spectrum of the magnetic field fluctuation in 
order to use Eq. El 

5. Simulated cluster radio sources depolarization 

In this section wc investigate the effects of cluster mag- 
netic fields on the depolarization of cluster radio sources. 

5.1. Beam depolarization of cluster radio galaxies 

This kind of depolarization is expected to be particularly 
strong if differential Faraday rotation is occurring within 
the beam, i.e. if the minimum magnetic field scale length 
is smaller than the beam size. 

Due to the high RM values in the inner regions of 
the cluster, the radio sources located in projection near 
the cluster center are expected to have a higher depo- 
larization than those located at the periphery. In what 
follows we refer to the depolarization, DP, as the ratio be- 
tween the polarization obtained after the rotation of the 
plane of polarization and the intrinsic polarization at a 
given frequency. To study the dependence of this effect 
on the cluster magnetic field strength and structure, we 
simulated the depolarization of virtual radio sources lo- 
cated halfway through the cluster at increasing distances 
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Fig. 4. Comparison of the ctrm profiles obtained from 
our simulations (solid line) and the analytical approxima- 
tion (dashed and dotted lines). We assumed a source sit- 
uated halfway through the cluster and the same magnetic 
field strength and gas density radial profiles in both the 
simulations and the analytical formulation. In the simula- 
tions we assumed a magnetic field power spectrum with a 
slope n — 2 and a range of scales from Amin = 6 kpc to 
Amax = 768 kpc. 



from the cluster center. The virtual radio sources are sup- 
posed to lie at the distance of the cluster center and their 
three-dimensional structure and inclination are neglected. 
Their intrinsic polarization intensity, Pint , was considered 
spatially uniform and completely ordered, i.e. the intrin- 
sic polarization plane direction, ^int, did not vary across 
sources. 

We first calculated the expected polarization plane di- 
rection, ^obs('^), at a given frequency along each line-of- 
sight of the sources according to Eq. ^ 

We then converted Pint and *obs into the Stokes pa- 
rameters Qiv) and U{v), by inverting the relations: 



P 



i arctan([//(3) 



(18) 



Finally, the simulated Q{v) and U{v) images have been 
convolved with the desired beam and transformed back to 
Pohs{v) and *obs(i^) using Eqs.CHl 

The parameters of the simulated cluster are those 
listed in Tab.^ In Fig.[Slwe show the depolarization of the 
sources at 1.4 GHz as a function of the distance from the 
cluster center for two different decreasing (/i=0.3) mag- 
netic field strength profiles with (B)o = 1 (top panels) 
and 10 jiG (bottom panels). 
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Fig. 5. Simulated 1.4 GHz beam-depolarization as a function of the distance from the cluster center in the case of 
a varying magnetic field (/z = 0.3) with a central strength of (B)o ~ 1 (top panels) and 10 /xG (bottom panels). The 
simulated cluster is placed at a distance of about 300 Mpc (z=0.05). Solid and dashed lines refer to a beam size of 
45" and 15", respectively. At the distance of the cluster these beam widths correspond to projected linear sizes of 
about 60 and 20 kpc. 



The simulated cluster is placed at a distance of about 
300 Mpc {z = 0.05). Solid and dashed Unes refer to a beam 
size of 45" and 15", respectively. At the distance of the 
cluster these beam widths correspond to projected linear 
sizes, Abeam, of about 60 and 20 kpc, i.e. Abcam/Amin — 10 
and Abeam/Ainin ~ 3, respectively. 

We found that with a beam size of 45" the dependence 
of the simulated polarization profiles on the adopted power 
spectrum spectral index is only marginal, while it becomes 
significant by considering a beam of 15". 

From Fig. it is evident that the depolarization ef- 
fect is strongly dependent on the cluster magnetic field 
strength. In the case of a weak magnetic field strength 
(e.g. 1 /iG) the Faraday rotation only affects the polariza- 
tion of the sources located in projection near the cluster 
center. As the average cluster magnetic field strength in- 
creases, the beam depolarization becomes significant at 
progressively larger distances. For a cluster magnetic field 
with a strength of (B)o — 10 fiG in the center, the sources 
are severely depolarized for distances up to 2-^3 re at both 
45" and 15" resolutions. 

The comparison of these simulations with the polar- 
ization properties of large samples of radio sources could 
provide an additional statistical constraint on the intra- 
cluster magnetic field strength. 

5.2. The morphology and (de)polanzation of the radio 
halo emission 

Cluster radio halos may provide important information 
about the spatial power spectrum of the magnetic field 
fluctuations on large scales (> 100 kpc). In particular, 
different values of the power spectrum spectral index will 
generate very different total intensity and polarization 
brightness distributions for the radio halo emission. 



For a given spectral index n, we calculated the ex- 
pected total intensity and polarization synchrotron emis- 
sion by introducing in the three-dimensional magnetic 
field an isotropic population of relativistic electrons whose 
distribution follows: 



N{e,9) ^Noe-%sine)/2 



(19) 



where e and 6 are the electron's energy and the pitch 
angle between the electron's velocity and the direction of 
the magnetic field, respectively. 

At each point of the computation grid, we calculated 
the synchrotron emissivity, Jsyn, by convolving the parti- 
cle energy distribution given in Eq. 1191 with the emission 
spectrum of the single electron 



de 
dtdv 



= CiB sin 9 F 



aSsin6le2 



where F{x) is the synchrotron radiation kernel 



oo 

F{x) = X j K^/'i{z)dz 



(20) 



(21) 



where C/ and C^, are constants and whose numerical val- 
ues (cgs units) are 6.26 x 10^^ and 2.34 x 10^^^, respec- 
tively, while ^^5/3 is the 5/3-order modified Bessel func- 
tion. 

Let i?_L be the magnetic field perpendicular to the line- 
of-sight. Because of the high beaming of the synchrotron 
radiation, only those electrons with = sm~^{B^/ B) are 
visible by the observer. Therefore we have 



^syn('^) — C*) 



N^B^ 
2 B 



F 



C,Bi_e 



e-^de 



(22) 



where 



and 



5 are respectively the low and high 
energy cut-offs of the electron distribution. We calculated 
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the intrinsic linear polarization emissivity, Pint , according 
to 



-Pmt(j^) = Jsynil^) 



I 

emin 


>fe 


V 


-y 


-'de 


Cm a- 


^g(^ 


y 

D 


-.1.- 


-Sde 



(23) 



where G{x) = xK2/3{x) (see e.g. Rybicki & Lightman 
1979). 

At our selected frequencies the radio halo is supposed 
to be optically thin, therefore we calculated the total in- 
tensity brightness distribution, /syn(i^), by simply inte- 
grating the emissivity in Eq. |211 along the line-of-sight 



'syn 



^--^I- 



Jsyn{v)dl 



(24) 



where L is the computational grid size. 

The calculation of the polarization brightness distri- 
bution must consider two important effects: i) since the 
direction of magnetic field changes from point to point 
inside the cluster also the orientation of the intrinsic po- 
larization plane is a spatially varying quantity and ii) as 
the radio waves propagate through the magnetized intra- 
cluster medium their polarization plane is subject to the 
internal Faraday rotation. The well known result of this 
effect is a frequency dependent depolarization of the in- 
trinsic synchrotron emission. 

We first calculated the rotated polarization plane di- 
rection, \l/(i/, /), of the emission from a volume element 
located at a position I along the line-of-sight (Z=0 corre- 
spond to the farthest side of the computational grid) : 



*(i.,0 = *i„t(0 + (cH'xRM(0 



(25) 



where ^int(0 i^ perpendicular to B±{1) and RM(Z) is 
the internal Faraday rotation occurring in between the 
volume element and the observer: 



RM(0 



812 / noB||d? . 



(26) 



It is then convenient again to convert Pint(i^) and 
'^{I'J) given in Eas. 051 and [^ into the Stokes parame- 
ters Q{v,l) and U{i',l), by inverting the relations given 
by Eqs.d 

Finally, wc obtained the simulated images of Q(i') and 
U{v) by adding together the contributes of all volumes 
along a given line-of-sight 

Q{iy) ^jQ{i^J)dl (27) 


L 

U{v) ^jU{v,l)dl 



These simulated polarization images take into account 
of both the effects of the field ordering and internal 
Faraday rotation. 



Determining the origin and evolution of the relativis- 
tic electron population responsible for the halo emission in 
cluster of galaxy is beyond the purposes of this work. Here 
we restrict our interest to the study the dependence of the 
radio halo morphology and polarization on the magnetic 
field strength and structure. We adopted the simplest rela- 
tivistic electrons distribution as possible which gives plau- 
sible values for the halo emission spectral index and total 
luminosity. 

In the simulation we adopt an electron energy spectral 
index (5 = 3. We choose a Lorentz factor corresponding to 
the high-energy cutoff of the electron distribution, 7max = 
emax/('7icC^), in a way that the high-frequency cutoff in 
the emission spectrum (at the cluster center) occurs at a 
frequency of 5 GHz: 



7n 



3.4 X 10"* {B) 



-0.5 



0.5 



iiG ''5GHz 



(28) 



The low-energy cutoff, 7min, and the relativistic electron 
density, A'o, are adjusted to guarantee a typical aver- 
age halo brightness at 1.4 GHz of /1.4 ~ 5 x 10^'* W 
Hz~^Mpc~^ (e.g. Giovannini et al. 1999) and a perfect 
equipartition between the relativistic electron and the 
magnetic field energy densities: 



7,„in ^ 1.9 X IQHB)-1^. 



(29) 



The relativistic electron density A^o follows the mag- 
netic field energy density decrease so that the equiparti- 
tion is satisfied at each distance from the cluster center. 

We simulated the expected total intensity and polar- 
ization brightness distribution at 1.4 GHz and 327 MHz, 
for three values of the magnetic field strength, namely 
(S)o =0.5, 1 and 5 ^G. 

Due to the assumed declining (/i = 0.3) mean magnetic 
field strength, the spectral index systematically increases 
from af ™H^ ~ 1 up to af/cif/ =^1.3 moving from the 
cluster center outward. 

So far, polarization emission from radio halos has not 
been detected. The current upper limits to the polariza- 
tion at 1.4 GHz are a few percent (3 — 4%) for beams of 
about 45". Due to the low spatial resolution of these obser- 
vations, the beam depolarization effect is expected to be 
very strong. Therefore, we simulated the polarization in- 
tensity expected at 45" resolution for a radio halo located 
at a redshift z — 0.05. At this distance, the correspond- 
ing linear resolution is of about 60 kpc. In the left-column 
panels of Fig. we show the simulated radio halo bright- 
ness and polarization percentage distributions at 1.4 GHz 
for three different values of the power spectrum spectral 
index n. In the right-column panels of Fig. we show the 
expected fractional polarization profiles at 1.4 GHz and 
327 MHz for the different values of the average magnetic 
field strength. 

As one would expect, we found that the degree of po- 
larization percentage: 

i) increases with increasing values of n; 
ii) increases with the frequency v] 
iii) decreases towards the cluster center; 



12 



M. Murgia ot al.: Magnetic Fields and Faraday Rotation in Clusters of Galaxies. 



fin «M W htr"-1 tJ^x'S^ 




■3 



ti^i 






" 




i.n dHi 




' 


- — 


.. E»r iHi> 








<n. 






r\ 


jiji f^ 




f\A 




,1 jiE 


v 


^ 




ptiid 


65s 




-_->i3 


;M-l-**i 



Dm -qm zm d -zm -qm -am 



h« H« iW MO 11« 



fiQ "m w i-te*-i wtc"^3 




£ 



n-a 


1 ■ ■ 1 


■ ' 


- 


. iA tlHl 


^ 





a^uHi 


<n>. 






^"■^ 




f 


W,.» 


v; 


'/ 


^,~,,-^ ^"^ 


T»Trw?r 


tf^^riiii 


-■^-^T H-»- — T 1 . . . 1 



tn 4U. «> 

B-isLetico (kpc-) 





^ 


1 




ID ai 


31 
1 1 


BID 
■VD 




OlD ■ 




J 


h 




p n. 

E 

-aiD - 




1 






^1 


-41D 
-BID ■ 


^^prr . . L-i^^j bis^. 


^j^H 


^^^^^^^^^^^^^^^^v^ 


J 



om ^m 211] D -zm -qtn -am 




tV« 4fC «* «0 LMfe 

l>lsLBnco (kpo) 



Fig. 6. simulated synchrotron halo emission for a cluster at a distance of z = 0.05 as it would be observed with a 
beam of 45" (~ 60 kpc). Left panels: simulated halo brightness and polarization images at 1.4 GHz for different values 
of the magnetic field power spectrum slope n and (-B)o =1 ^G; the vectors lengths are proportional to the degree 
of polarization, with 100 percent corresponding to 100 kpc on the sky. Field directions are those of the E-vector. 
Right panels: radially averaged profiles of the polarization percentage at 327 MHz and 1.4 GHz for three values of the 
magnetic field strength, namely (i?)o =0.5, 1 and 5 /iG. 
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iv) decreases with increasing magnetic field strength. 

Our results indicate that a power spectrum slope 
steeper than n — 3 and a magnetic field strength lower 
than ^ 1/iG result in a radio halo polarization percent- 
age at a frequency of 1.4 GHz that is far in excess the 
current observational upper limits at 45" resolution. This 
means that either the power spectrum spectral index is 
flatter than n = 3 or the magnetic field strength is signifi- 
cantly higher than ~ 1/iG. The halo depolarization at 327 
MHz is particularly severe and the expected polarization 
percentage at this frequency is always below 1%. 

We also found that the magnetic field power spectrum 
slope has a significant effect in shaping the radio halo. 
In particular, flat power spectrum indexes {n < 3) give 
raise to very smooth radio brightness images (under the 
assumption that the radiating electrons are uniformly dis- 
tributed). RM images of radio galaxies in conjunction with 
observations of the halo brightness and polarization dis- 
tributions could provide information on the behavior of 
the non-thermal (magnetic field and relativistic electrons) 
components of the intra-cluster medium over a large range 
of spatial scales. For this reason, it is important to apply 
our approach to those cluster with both types of radio 
sources, such as A2255, in order to determine if a single 
magnetic field power spectrum can account for both the 
observed RM and radio halo images (Govoni et al. 2002, 
Govoni et al. 2004 in preparation). 



6. Application to the data 

6.1. The case of Abell 119 

The cluster of galaxy A119 is an ideal study case to apply 
our new approach. It is characterized by three extended ra- 
dio galaxies, 0053-015, 0053-016 and 3C29, located at dif- 
ferent projected distances from the cluster center, respec- 
tively at 170, 453 and 1515 kpc (see Tab. Eland Fig.HTIkl. 
All these sources are highly polarized and they have been 
studied in detail with the Very Large Array (VLA) by 
Feretti et al. (1999a). From RM images with a resolution 
of 3.75"(~ 4.4 kpc) Feretti et al. (1999a) estimated a tan- 
gling scale of ~ 10 kpc and, by using the simple analytical 
approach (see Sect. l4.4|l . they inferred a strength of about 
5/iG for an assumed constant intra-cluster magnetic field. 
However, Dolag et al. (2001) by comparing the X-ray sur- 
face brightness with the ctrm of the three radio galaxies 
(see Sect. 3.1) estimated that the magnetic field strength 
in A119 scales as i? oc n°-^. Taking into account such a 
varying magnetic field, the analytical approach yields a 
central magnetic field of ~ 7.5 /iG. 

We used FARADAY in order to numerically obtain the 
expected RM images for the radio galaxies of Abell 119. 
In the numerical simulations we used the gas density dis- 
tribution deduced by X-ray observations for this cluster 
(re = 378 kpc, tt-q = 1-18 x 10^^ cm'^ and /? = 0.56; 
Cirimele et al. 1997), and we left the magnetic field power 
spectrum slope and strength as free parameters. All three 



radio galaxies are assumed to lie on a plane perpendic- 
ular to the line-of-sight at the distance of cluster center. 
The details of the three-dimensional structure of the ra- 
dio source are neglected, and the entire observed Faraday 
rotation is assumed to occur in the cluster intergalactic 
medium. 

A119 is located in a region of low Galactic RM (see 
Tab. |2Jl and therefore we apply no correction to the ob- 
served (RM) with this respect. We compared both ctrm and 
(RM)| versus the X-ray surface brightness S^ taken from 
Dolag et al. (2001) to estimate the index fi for the radial 
scaling of the magnetic field. Fig. [3 shows the same de- 
pendence RM oc S^'^ for both quantities, in agreement 
with the result by Dolag et al. (2001). From Eq.[l3| with 
/3 = 0.56 and a = 1.1, it follows that /i ~ 0.5. 

As a first step we compared the observed ctrm and 
(RM)| with the expectation of the simulations. Since 
the linear resolutions of both the computational grid and 
the observations are very similar, we do not apply any 
Gaussian smoothing to the simulated RM images. Fig. |21 
show the best fit of the simulated ctrm radial profiles for 
the radio galaxies in A119. The typical uncertainties in 
the observed RMs are of the order of 10—15 rad m~^. 
However, a realistic error budget should also consider the 
effects of a number of systematic uncertainties which are 
difficult to estimate, such as the relative distance of the 
galaxies from the cluster center along the line-of-sight 
and the statistical fluctuations which are intrinsic to the 
random nature of the magnetic field. Moreover, the "ob- 
served" RMs themselves are the result of a fit procedure 
as outlined in Sect. 2. For all these reasons, we decided 
to use a simplified approach, i.e., a non- weighted least 
square fit without attempting to include either ctrm nor 
(RM) uncertainties in the fit procedure. 

In the fit procedure we fixed n at a given value and let 
the magnetic field strength as a free parameter. The best 
fit of the magnetic field strength, obtained by fitting the 
crm profiles (Fig. |S1 left), gives (B)o — 5 /xG for each n. 
This is the mean magnetic field at the cluster center. Since 
the modeled magnetic field strength decreases from the 
cluster center outward (^=0.5), the mean field strength 
within a spherical volume of radius r = 3rc (with r^ — 
378 kpc) gives (B)~ 1.5 fiG. The comparison between 
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Fig. 8. Best fit of the simulated RM ctrm profiles for the cluster A119 (left column panels). The comparison of 
simulated and observed |(RM}| is shown in right column panels. Solid dots refer to the observed values for the radio 
galaxies 0053-015, 0053-016 and 3C29. The best fit to the central magnetic field strength gives (B)o — 5 /iG. Since the 
modeled magnetic field strength decreases from the cluster center outward (/z ~ 0.5), the mean field strength within a 
spherical volume with radius r = 3 re is (B)~ 1.5 iiG. We find that the spectral index n = 2 gives the best description 
of the data if the magnetic field power spectrum ranges from a minimum scale of 6 kpc up to a maximum scale of 768 
kpc. The fit corresponding to the models with n < 1 and n > 2 are definitely poorer since to the |(RM)|/(Tr,m ratios 
they predict are too low and too high, respectively. The Kolmogorov spectral index (n = 11/3) gives a reasonable fit 
to the data only if the maximum scale of the magnetic field fluctuations is about 100 kpc (bottom panels). 



simulated and observed |(RM)| (Fig. |H1 right) allows us 
to constrain the value of the index of the power spectrum 
to n ~ 2. Fits corresponding to models with n > 2 are 
definitely poorer as evident from the high |(RM)| values 
they predict. Also models characterized by n < 2 are a 
poorer fit because their low expected |(RM)| values are 
incompatible with the observations. As shown in bottom 
panel of Fig.|Sl the Kolmogorov spectral index (n = 11/3) 



gives a reasonable fit to the data only if the maximum 
scale of the magnetic field fluctuations is about 100 kpc. 

As expected from the discussion in Sect 14. 31 the ob- 
served ijRM and the |(RM)| profiles, can be explained ei- 
ther by a broad and flat {n ~ 2) or by a narrower but 
steeper {n ~ 11/3, Amax^lOO kpc) magnetic fleld power 
spectrum. 
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Fig. 9. Comparison of observed and simulated RM images of the source 0053-016 in A119. The observed RM image 
is shown in the top left panel. Models with n ~ 2 better reproduce the observed RM features. This is confirmed by 
the structure functions of the observed and simulated RM images shown in the bottom left panel. 



A way to discriminate the two models is to compare 
the structure functions of the observed and simulated RM 
images of the individual radio galaxies. As an example, in 
Fig. 121 we show the observed and simulated RM images 
relative to the source nearest to cluster center: 0053-015. 
The model with n = 2 is in remarkable agreement with the 
observed RM features. This is confirmed by the structure 
functions'* plotted in the bottom left panel of Fig. |51 The 
observed RM structure function cuts off below a linear 
scale of about 10 kpc. This cutoff and the overall shape of 



* To calculate the structure function of the RM images we 
first construct the structure function RM matrix 



S{dx, dy) = ([RM(x, y) - RM(2: + dx,y + dy)Y) 



(30) 



where dx and dy are both positive and negative offsets and the 
ensemble average is calculated over all the positions, {x,y), in 
the RM image. The RM structure function is then obtained 
by radially averaging S{dx, dy) in annuli of increasing size A — 
^/dx^ + dy^. 



the observed RM structure functions is reproduced well by 
the model with n = 2. We note that even if the magnetic 
field power spectrum with the Kolmogorov index n — 11/3 
and Ainax = 100 fits the observed RM profiles, it generates 
a structure function which is not in agreement with the 
data. We found similar results for the other two sources 
(not shown). 

6.1.1. The composite RM image 

To confirm the above results we calculated the two- 
dimensional power spectrum of a composite RM image of 
A119 obtained by combining all together the data of the 
three radio galaxies of the cluster (Fig llOb ). In Fig. llOb 
we shown the power spectrum of the "mask" , i.e. a win- 
dow function image which is equal to unity if the line-of- 
sight intercepts one of the sources and zero otherwise. The 
power spectrum of the composite RM image is shown in 
Fig. llOb . This is the convolution of the true power spec- 
trum with the power spectrum of the mask. Since the 



16 



M. Murgia et al.: Magnetic Fields and Faraday Rotation in Clusters of Galaxies. 



a) 



Qja upc. 



iJSHpc 



^^^^^^^^^^~^ UHpc 





d) 



W-^ ^ 



IRMJ* 



10" - 



lfjl« - 



10*! . 



ID 



k (1/«pc) 

: — . 1 I Mil 



ICI^^ 



DlncniBd 



. ^^^_ ii=S. 



A„,-]PO S(pc 




Al 13 



I i • 



lOflO 



li I J-MpE| 



1DQ 10 

A(kpc) 



I L 



Fig. 10. a) Composite RM image of A119 obtained by combining together the individual images of the three radio 
galaxies. The cross marks the position of the X-ray centroid. b) The power spectrum of the mask window function. 
The short-wavenumber modes are due to the interference between the source-to-source baselines. The inset shows the 
radially averaged profile of the mask power spectrum, c) Power spectrum of the composite RM image, d) Comparison 
of the observed (solid line) and simulated (dashed lines) radially averaged RM power spectra, see text. 



power spectrum of the mask is quite broad, we cannot 
measure the power spectrum slope directly from the im- 
age shown in Fig. llOb . However, we can multiply by the 
same mask the simulated RM image corresponding to the 



models with n = 2 and n — 11/3 with Amax = 100, calcu- 
late their power spectra, and compare them with the ob- 
served one. The radially averaged profiles of the observed 
(solid line) and simulated (dashed lines) power spectra of 
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the composite RM image are shown in FigllUli. Again, we 
find that the model with n = 2 gives the best description 
of the data. 

To summarize, the best fit of both the RM profiles, 
structure functions and power spectra observed for the 
three extended radio galaxies in A119 is obtained for a 
spectral index n — 2, Ai„in=6 kpc, Ainax=770 kpc and 
(B)o — 5 /xG. The mean field strength within a spheri- 
cal volume of radius r = 3rc (with re = 378 kpc) gives 
(B)~ 1.5 ^G. The observed |(RM)|/o-rm ratios rule out 
n ~ I while models with n > 2, compatible with the 
constraint A^ax — 100 kpc, produce RM images whose 
structure functions and power spectra are incompatible 
with the data. These results are based on the assumption 
that the observed Faraday rotation is produced entirely 
by a magnetized intra-cluster medium. 

6.2. Variations of the power spectrum spectral index 
among clusters 

We have seen in the case of A119 how by comparing the 
simulated RMs with the observed RMs obtained for radio 
sources located at different projected distances from the 
cluster center, it is possible to gain important information 
about the magnetic field power spectrum and strength. 
However, one important question concerns the variation 
of the magnetic field power spectrum among clusters of 
galaxies. To explore this issue, we collected a sample of 
nine clusters and groups of galaxies (without cooling flows) 
for a total of 14 radio sources for which high-quality RM 
data are available in the literature. The cluster sample we 
chose for this study is listed in Tab. 13 Since we are inter- 
ested in the mean and dispersion of the RM, we included 
in this analysis only those clusters hosting radio sources 
for which high-resolution RM images are available. 

The clusters in our sample are characterized by quite 
different gas density profiles, linear size and probably mag- 
netic field strengths. However, since these parameters af- 
fect both |(RM)| and ctrm in the same way, it is possible 
to combine the RM data from different clusters and use 
the ratio |(RM)|/crRM to study the variation of the spec- 
tral index n among them from a statistical point of view. 

The mean of the RM is a fundamental parameter for 
our analysis and therefore the contribution of our own 
Galaxy to this quantity should be removed carefully. We 
determined the galactic contribution to the (RM) for each 
cluster of our sample according to the values reported 
by Simard-Normandin et al. (1981). We use the distance- 
weighted mean of the RMs of all sources located within 
±15deg of galactic longitude and latitude, respectively. 
The galactic RM contributions are listed in Tab. [21 (col. 
4). 

Fig. E] shows the trend of the observed 
|(RM)|/(Tr{M ratios as a function of the sampling re- 
gion size (As) for n = 1, 2, 3 and 4. The observed values 
of the |(RM)|/(Trm ratio, lie above 0.1 and below 2 for 
all clusters except for 3C 449, which has a very low 
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Fig. 11. The ratio |(RM)|/ctrm as a function of the sam- 
pling region. The horizontal lines represent the simulated 
trends corresponding to n—1, 2, 3 and 4. 



(RM)| if corrected for the galactic RM contribution. 
It is important to note that 3C 449 belongs to a small 
group and is located at low galactic latitude, therefore 
the correction for the galactic RM contribution should be 
considered carefully since it could strongly influence the 
result. The mean of |(RM)|/crFiM ratio calculated for the 
14 radio sources is ~ 0.6. 

The values of |(RM)|/crRM expected from the simu- 
lations for n > 3 exceed those observed by a significant 
amount. Indeed, as already found for A119, the rather 
flat spectral index n ~ 1 — 2 represents the best descrip- 
tion of the data. The relatively large dispersion of the ob- 
served data is largely consistent with what expected from 
the intrinsic statistical scatter shown in the right panel of 
FigEl (the scatter is larger towards the lower values of the 
(RM)|/ctrm ratio). 

In A2255, for both radio galaxies the mean of the ob- 
served RM is greater than its dispersion: |(RM)|/crFiM— 
1.2. This is a strong indication that, at least for this clus- 
ter, the magnetic field power spectrum might be steeper 
than in the other cases. The inferred power spectrum index 
could be as high as n = 2 or more. A2255 is a very interest- 
ing cluster since it also contains a wide synchrotron radio 
halo and the magnetic field power spectrum slope has a 
major role in determining the large scale morphology of 
such a diffuse emission (see Sect. 15. 2f) . 
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Table 2. RM data taken from the literature. 



Cluster 


source 


Distance 
(kpc) 


Gal.RM 
(rad/m^) 


(RM) 
(rad/m^) 


0"RM 


Ref. 


Coma 


NGC4869 


192 


-6 


-127 


181 


1 


A119 


0053-015 

0053-016 

3C29 


170 
453 
1515 


-1 


-F28 
-79 

+4 


152 
91 
13 


2 


3C129 


3C129.1 
3C129 


50 
500 


-10 


-^21 
-125 


200 
82 


3 


A514 


PKS0405-203 
PKS0446-206 


2098 
750 


-3 


-19 

46 


48 
47 


4 


A2255 


B1712+640 
B1713+641 


392 
612 


-6 


-91 

+67 


78 
59 


5 


Zw0056.9+2636 


NGC326 E-lobe 
NGC326 W-lobe 


70 
35 


-36 


-16 
-21 


58 
60 


6 


A400 


3C75 





4 


-7.6 


100 


7 


A2634 


3C465 





-60 


-25 


120 


7 


3C449 


3C449 





-165 


-166 


38 


8 



Column 1 : Cluster name; 2: Source name; 3: Source projected distance from cluster center; 4: Galactic RM contribute; 5: RM 
mean; 6: RM dispersion; 7: References; 

References: 1) Feretti et al. 1995; 2) Feretti et al. 1999a; 3) Taylor et al. 2001; 4) Govoni ct al. 2001a; 5) Govoni et al. 2003 in 
prep; 6) Murgia et al. 2001; 7) Eilek & Owen 2002; 8) Feretti et al. 1999b. 



7. Summary 

In this paper we present a numerical approach to investi- 
gate the strength and structure of magnetic fields perme- 
ating the intergalactic medium in clusters of galaxies. The 
estimates of cluster magnetic field strengths appearing in 
the literature have been calculated through their effects 
on the polarization properties of radio galaxies in or be- 
hind clusters by using an analytical formulation based on 
the approximation that the magnetic field is tangled on a 
single scale. However, detailed observations of these radio 
sources and MHD simulations suggest that it is necessary 
to consider cluster magnetic fields which fiuctuate over a 
wide range of spatial scales. To accomplish this we sim- 
ulated a random three-dimensional magnetic field with a 
power law power spectrum in a cubical box. We calcu- 
lated the effects of the expected Faraday rotation on the 
polarization properties of radio cluster sources such as ra- 
dio galaxies and radio halos whose radiation crosses the 
magnetized intracluster medium. The main issues of this 
paper arc summarized as follows: 

• We determined the statistical behavior of the RM 
over a typical radio galaxy-sized region as a function of 
distance from the cluster center for different values of the 
magnetic field power spectrum slope and range of scales. 
The |(RM)|/crr{M ratio is an easily measurable observable 
that is expected to refiect the relative amplitude of large 
to small-scale magnetic field fluctuations. The numerical 
simulations allow a quantitative evaluation of this effect. 
We found that the |(RM)|/crFiM ratio depends rather criti- 
cally on the power spectrum slope, n, and therefore it can 
be used to constrain this parameter if spatially-resolved 



RM images of radio sources located at different distances 
from the cluster center are available. 

• The magnetic field strength estimated by our sim- 
ulations in the case of n = 2 (i.e. most of the magnetic 
field energy resides in the small scales) results in a factor 
of about 2 lower strengths than the expectations based 
on the single-scale magnetic field approximation (see Fig. 
4). The analytic formulation leads to a reliable estimate 
for the average magnetic field strength only provided that 
the magnetic field correlation length is adopted (i.e. the 
magnetic field power spectrum is known). 

• We simulated the expected dependence of the sta- 
tistical beam depolarization of marginally resolved radio 
sources located at increasing distances from the cluster 
center. Wc found that the amount of depolarization is 
strongly dependent on the magnetic field strength. 

• We point out that radio halos may provide important 
information about the spatial power spectrum of the mag- 
netic field fiuctuations on large scales. In particular, differ- 
ent values of the power spectrum spectral index give rise 
to very different total intensity and polarization bright- 
ness distributions. So far, polarization emission from ra- 
dio halos have not been detected. Our results indicate that 
a power spectrum slope steeper than n = 3 and a mag- 
netic field strength lower than ^ IfiG predict a radio halo 
polarization percentage at a frequency of 1.4 GHz that 
is far in excess the current observational upper limits at 
45" resolution. This means that either the power spec- 
trum spectral index is fiattcr than n = 3 or the magnetic 
field strength is significantly higher than r^ 1/iG or both. 
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• In the last part of the paper, simulations and data 
have been compared through an interactive approach. The 
cluster A119 is an ideal case for our new approach. It is 
characterized by three extended radio galaxies located at 
different projected distances from the cluster center and 
it is located in a region of very low Galactic RM. The best 
fit of both the RM profiles and image structure functions 
observed for the three extended radio galaxies in A119 is 
obtained for a spectral index n = 2 and (B)o — 5 /iG. 
This is the mean magnetic field at the cluster center, the 
mean field strength within a spherical volume with radius 
r = 3rc is (B)~ 1.5 /iG. The observed |(RM)|/crRM ratios 
rule out n — 1 while models with n > 2, compatible with 
the constraint A,„ax — 100 kpc, produce RM images whose 
power spectra are incompatible with the data. 

• An as-yet unexplored issue concerns the variation of 
the magnetic field power spectrum among cluster of galax- 
ies. For this purpose, we collected a sample of nine clus- 
ters (without cooling flows) for a total of 14 radio sources 
for which high-quality RM data are available in litera- 
ture. The values of |(RM)|/crp{M expected from the simu- 
lations for n > 3 exceed the observed ones by a significant 
amount. Indeed, as already found for A119, the rather flat 
spectral index n ~ 2 represents the best description of the 
data. 
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